Zero Sound Propagation in Femto-Scale Quantum Liquids 
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Charge equihbration has been recognized as a dominant process at the early stage of low-energy 
heavy-ion reactions. The production of exotic nuclei is suppressed under the appearance of charge 
equilibration, in which the proton-neutron ratios of the final reaction products are inevitably aver- 
aged. Therefore charge equilibration plays one of the most crucial roles in the synthesis of chemical 
elements. Focusing on how and when the charge equilibration takes place, the zero-sound propaga- 
tion in femto-scale quantum liquids is explained. 
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INTRODUCTION 

This review article is concerned with Fermionic prop- 
erty of heavy ions (many-nucleon systems with the size 
up to several 10-femtometers) colliding at the energy of 
a few MeV per nucleon. In the following, we refer to 
those collisions as low-energy heavy- ion collisions. Fusion 
is not necessarily achieved in the low-energy heavy-ion 
collisions, neither is fragmentation. Reaction dynamics 
and the resulting products can be drastically different de- 
pending on the impact parameter, the mass of colliding 
ions, and so on. Therefore, the reaction mechanism of 
low-energy heavy-ion collisions is worth investigating to 
understand the production of chemical elements. This is 
deeply concerned with an open problem as for the exis- 
tence and the origin of chemical elements including their 
production mechanism. 

Charge equilibrium in heavy-ion reactions means the 
states with the proton-neutron ratio corresponding to the 
average of the two colHding ions, and charge equilibration 
is the process leading towards charge equilibrium. The 
chemical property of final products can be different de- 
pending on whether charge equilibration appears or not, 
being determined by the number of protons. There is a 
relatively long research history for the charge equilibra- 
tion. In fact many experiments of low-energy heavy-ion 
collisions were carried out in the 1960's. In those ex- 
periments most of the final products were reported to 
be in charge equilibrium, even when the fragmentation 
takes place There should not be any restrictions 

to the proton-neutron ratio of the final products if the 
final product forms a kind of stable bound system, so 
that these experimental results cannot be trivially un- 
derstood. 

One of the most important features of charge equilibra- 
tion is its rapidness taking as much as a few 10~^^ s [ij. 
This time scale is actually short in the order of mag- 
nitude compared with the typical reaction time of low- 
energy heavy- ion collisions (1000 fm/c ~ 10~^^ s), so 
that charge equilibration has been recognized as an in- 
evitable and dominant process in low-energy heavy-ion 
reactions. The relation between charge equilibration and 



the isovector giant dipole resonance has been studied 
relatively well because of the correspondence in their 
time scales, but no decisive conclusion has ever been ob- 
tained. Indeed, including the question of "when does 
charge equilibration take place ?" , many things could not 
be explained merely by the isovector giant dipole reso- 
nance. For the theoretical research on the relation be- 
tween charge equilibration and the isovector giant dipole 
resonance, the importance of dipole mode to charge dis- 
tribution of a fissioning nucleus was pointed out by Hill- 
Wheeler in the 1950's. Research on charge equilibra- 
tion using time-dependent mean field calculations were 
started in the 1970's (Bonche-Ngo [3]). Recently re- 
search based on three-dimensional time-dependent mean 
field calculations was carried out by Simenel-Chomaz-de 
France 0, HI- In this article, based on Ref. [HI], the un- 
known relation between charge equilibration and the zero 
sound propagation is presented. It leads to a rather uni- 
versal recognition of zero sound propagation in femto- 
scale quantum liquids, which cannot necessarily be re- 
duced to the giant dipole resonance. 



THEORY OF MATTER WAVE PROPAGATION 
Landau's Fermi liquid theory 

Let us denote a many-particle wave function by 
^(t, xi, • • • , Xn). It is assumed to satisfy 



ih 



dt 



where H denotes the Hamiltonian operator. The solu- 
tion can be represented by ^(t^x) = e~*^^/^^(0, x) un- 
der a suitable boundary condition, if H does not depend 
on ^Xn)' Let the corresponding probability 

density be p{t,x), and begin with the classic theory of 
sound propagation inside gases. Readers may wonder 
why sound propagation is related to charge equilibration, 
it will be clarified step by step. Let the equilibrium prob- 
ability density be po{x). If the fluctuation is added to the 
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equilibrium: 



a force arises from the gradient pressure. Here the essen- 
tial property of sound propagation is extracted from the 
simplified linearized analysis. The equation of motion is 
given by 



(1) 



where v and P denote the velocity and the pressure, re- 
spectively. The pressure P depends on both p and the 
entropy S. On the other hand, the equation of continuity 
is given by 



d8p 
dt 



-div{pv). 



The right hand side is approximated by — po div v to 
the lowest order. After differentiating this equation by 
an equality is derived together with the divergence of Eq. 
([1]) . If we further assume that P is expanded with respect 
to (5p at a fixed entropy, the lowest order contribution 
brings about 
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This is a wave equation for 6p^ where ^/ {dP/dp)^ repre- 
sents the propagation speed of the density. For example, 
^y {dP/dp)g is given by vp/V^ for the perfect Fermi gas 
in its ground state {S = 0), where vf means the Fermi 
velocity. This type of propagation is called first sound, 
which provides a picture for particles propagating with 
changing density. 

On the other hand, we should pay attention to the 
propagation of particle without changing density. This 
type of propagation is called zero sound, which has a 
finite frequency even when the wave number is equal 
to zero (cf. zero point vibration). To understand zero 
sound, the linear response is considered for a given ex- 
ternal field Hex{t): 

Hex{t) = j d^x n{t,x) Uex{t,x), 

where n(t, x) denotes the number density operator. In 
particular, if we restrict ourselves to the impulsive per- 
turbation Uex{t^x) = /7ea^e*^^(5(t), the linear response is 

6 < n{t^ x) > 



UJ 



where tv^{q,(ju) denotes the retarded generalized dielec- 
tric function. The resonance frequency can be calcu- 
lated by the pole of the integrand of the right hand 



side. Eventually we assume the phonon dispersion re- 
lation {^}q = Cog), which reproduces the frequency of 
propagating wave, and is nonzero even when the wave 
number is equal to zero: 

l = V{q)U°''{q,ng-ng), 

then the resonance frequency is obtained. Consequently 
the dispersion relation at long wave length becomes 
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1, 



where x = muo/hkFq] m and kF denote the mass and the 
Fermi wave number, respectively. Note that this relation 
is the representation in momentum space. The lower- 
limit of the velocity for zero sound is the Fermi velocity, 
because the non-damping mode can only exist when x > 
1 (see the denominator inside the logarithmic function). 
If V{q) V{0) is assumed in the limit g ^ 0, the left 
hand side becomes j (mi?V^(0)), and the propagation 
velocity of zero sound is represented using F(0). For the 
weak-coupling limit (1^(0) << h? /rnkF)'- 
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mkFV{0) 

and for the strong-coupling limit (V^(O) >> /mkF)'- 
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This shows that the propagation velocity of zero sound is 
almost equal to in the weak-couphng limit. It is just 
VS times faster than the previously seen first sound veloc- 
ity of the perfect Fermi gas. A rather general discussion 
shows that zero sound is faster than first sound [il, 0]- 
Roughly speaking, the matter waves without changing 
density propagate more easily than those with changing 
density, because the change in density possibly leads to 
the appearance of larger restoring force. Zero sound pro- 
vides a picture of particle exchange within a quite short 
time (but not instantaneous). 

Two different sound propagations in Fermi liquids have 
been discussed in terms of whether they entail the den- 
sity change or not. In particular there exists zero sound 
in Fermionic many-body system, which is different from 
and faster than the ordinary sound. Zero sound has 
been known to arise from the collective dynamics of the 
Fermionic many-body systems. All the collective dynam- 
ics is actually based on the propagation of zero sound, 
while there are various representations (various modes) 
for the collective dynamics. One distinct difference be- 
tween zero and first sounds is their relation to collisions 
between particles. First sound appears when the states 
are in local thermal equilibrium. This corresponds to 
the situation when the mean interparticle collision time 
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is sufficiently smaller than the oscillation period of the 
propagating wave. Meanwhile zero sound is associated 
with the collective excitation mode, which disappears 
when there are many collisions between particles. This 
corresponds to the situation when the mean interparti- 
cle collision time is larger than the oscillation period of 
the propagating wave. Although such a collisionless en- 
ergy regime is expected to be realized due to the Pauli 
principle, its validity should be conffimed for an indi- 
vidual physical system. Consequently zero sound is ex- 
pected to be important to low-energy phenomena, while 
first sound becomes effective at higher energies. For the 
usage of the terminology low and high energies, we had 
better note that both sounds are scaled hy vf (i.e. Fermi 
energy). For more details of sound propagation, see Sec. 
5 of Ref. S. 



Nucleon propagation in heavy-ion collisions 

As a theoretical research on the zero sound in 
Fermionic many-body systems. Landau's Fermi liquid 
theory [qI is well known, where zero sound propagation 
was actually seen in liquid ^He [lo|- However, there is 
no guarantee that such a sound plays a role in heavy-ion 
reactions. There are two essential differences between 
heavy- ion reactions and the ^He case; the physical sys- 
tem consists of finite numbers of nucleons, and the event 
is accomplished within a finite time interval. That is to 
say, both size and time are highly restricted in heavy-ion 
reactions. Indeed, in the context of many-nucleon sys- 
tems, the main interest of zero sound propagation was 
not in heavy-ion reactions but in nuclear vibrations (for 
example, see Ref. [11]). In the following the mechanism of 
charge equilibration is discussed with respect to whether 
it is achieved by nucleon propagation with or without 
changing density. 



MECHANISM OF CHARGE EQUILIBRATION 
Nucleon propagation realizing charge equilibration 

As an example of Fermionic many-body systems, we 
consider a many-nucleon system. There are two types of 
nucleons, that is, protons and neutrons. Charge equili- 
bration is the mixing of protons and neutrons due to the 
time evolution, therefore it is a kind of chemical mix- 
ing. Apart from such a chemical equilibration, there 
are several kinds of equilibration in heavy- ion reactions, 
that is, mass equilibration, momentum equilibration, and 
thermal equilibration. Among them charge equilibration 
(chemical equilibration) has drawn special attention be- 
cause of its crucial role in the synthesis of chemical ele- 
ments. 



Let us consider the collision of two ions. The existence 
of an upper energy limit for charge equilibration has been 
pointed out by Iwata-Otsuka-Maruhn-Itagaki [5] for the 
first time. The limit is presented by a formula explain- 
ing both experiments and numerical calculations based 
on microscopic three-dimensional time-dependent mean- 
field theory. Furthermore, the upper energy limit has 
been concluded to be determined by the Fermi energy. 
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where i = 1,2 identifies the two colliding ions, and 
Eiab means the energy in the laboratory frame. Ai and 
^2, which satisfy A = Ai ^ A2, Ai = Zi + TVi and 
A2 = Z2 -\- denote the masses of the two colliding 
ions, where Zi, Z2, A^i and A^2 denote the proton and 
the neutron numbers of each colliding ion (labeled by z), 
respectively, e and S are parameters introduced based on 
Ref. [l2| . This formula provides an upper energy limit of 
charge equilibration, which arises from the nucleons prop- 
agating at the Fermi velocity. In heavy- ion collisions, 
there are four different Fermi velocities (Fermi energies), 
because there are two different kinds of nucleons and two 
colliding ions. The first term on the right hand side of 
Eq. (j3j) calculates those four Fermi velocities. In par- 
ticular, for the many-nucleon system, special structures 
such as skin and halo exist, and those structures change 
the Fermi energy. These effects are taken into account in 
e and S in Eq. (|4]). For the nucleon wave associated with 
each Fermi energy to propagate throughout both collid- 
ing ions after touching, the minimum is taken in Eq. (j4]). 
This treatment does not seem to be so important at a 
glance, but it contributes to derive the fact that different 
reactions having exactly the same composite nucleus can 
result in a different upper energy limit (for the difference 
of upper energy limit, see Fig. 2 of Ref. [HI). If it is not 
for such a treatment, we cannot explain both numerical 
calculations and experiments. Furthermore, for a correct 
comparison, it is necessary to estimate the relative veloc- 
ity at contact time. The most considerable effect here is 
the deceleration due to the Coulomb repulsion, which be- 
comes prominent for cases when the masses of the collid- 
ing ions are larger. This effect is considered in the second 
term on the right hand side of Eq.(j3]). Consequently it 
was confirmed in Ref. [HI that the final fragments achiev- 
ing charge equilibrium drastically decrease at an incident 
energy higher than the energy shown in the formula. 

The nucleons propagating at the Fermi velocity, which 
is represented by Eqs. ^ and (j4]), correspond to the 
propagation of zero sound. Indeed, according to the 
calculation of the mean free path of nucleons inside nu- 
clear matter [l3| , collisions between nucleons were shown 
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FIG. 1: (Color online) Propagation of neutron-rich flow is depicted for the collision between ^^Ca and ^^Ca, where ^^Ca and 
^^Ca correspond to the ions coming from the left and right hand sides, respectively. The colored parts correspond to the parts 
with N/Z > 1 (each frame is 40x30 fm^), and the density contour equal to 0.02 fm~^ is shown by a thick black curve. Three- 
dimensional time-dependent Hartree-Fock calculations with a Skyrme interaction (SLy4d) is carried out; the single-particle 
wave functions are represented on a Cartesian grid with the spacing of 0.8 fm, and the time unit of calculation is set to 1.5 
xlO~^^ s. The initial distance between the two colliding ions are set to 20 fm, then the relative velocity of collision is given. 



to appear and increase rapidly if the incident energy of 
heavy- ion collisions becomes a few 10 % larger than the 
Fermi energy (of the composite nucleus). In such situa- 
tions, zero sound itself disappears. The possible propa- 
gation speed of zero sound is not not different from the 
amplitude of the Fermi velocity. This fact explains the 
reason why Eqs. (j3j) and (jH are related to zero sound 
propagation. 

Figure [1] visualizes the propagation of charge equili- 
brating flow: neutron-rich flow from ^^Ca (N/Z = 32/20) 
to ^^Ca (N/Z = 16/20). Note here that, because this 
heavy-ion reaction is classified to class II of the classi- 
fication shown in Fig. 1 of [1^ and the "N/Z= l"-line 
is located between the initial points of ^^Ca and ^^Ca 
on the N-Z plane, its charge equilibration dynamics can 
be measured only by the flow of neutrons. Cases with 
two different impact parameters and two different ener- 
gies are compared. The final product are different, where 
only the case with 6 = 6.1 fm and E/A =1.0 MeV results 
in fragmentation, while the other cases result in fusion. 
Despite such significant differences in final products, the 
propagation speed of charge equilibrating flow is almost 
the same. It can be confirmed by the time evolutions of 
0.75 X 10~^^ s after the neck formation (the neck is formed 



at 6.0x10-^^ s for E/A =1.0 MeV cases and 3.0x10"^^ s 
for E/A =2.0 MeV cases). It shows that the propagation 
speed of charge equilibrating flow is mostly independent 
of the incident energy (the relative velocity of collision) 
and the impact parameter, while the total contained neu- 
trons in neutron-rich flow depends highly on the incident 
energy. In particular the propagation speed of charge 
equilibrating flow is faster than the relative velocity of 
collision (Tabled]). 

With respect to zero sound, what was clarified in 
Ref. [il can be summarized in the following three points. 
First, zero sound propagation plays a role in heavy- ion 
collisions. Second, the fast charge equilibration, which is 
achieved within the order of 10 s, is realized as nu- 
cleon propagation without changing density. Third, there 
exists an upper energy limit for the fast charge equilibra- 
tion, which corresponds to the energy limit at which zero 
sound can survive. For the terminology of "fast" charge 
equilibration, it takes into account the existence of an- 
other kind of charge equilibration that has nothing to do 
with zero sound propagation. Such charge equilibration, 
which appears at higher energies, is more related to the 
first sound. Therefore its process is relatively slow com- 
pared to the fast charge equilibration, and insufficient 
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TABLE I: Comparison of speeds, where \vf\ is fixed to 1/3 of the speed of light (corresponding to the nuclear standard value). 
The propagation speed of charge-equilibrating flow is calculated by the propagation speed of the wave front of N/Z = 1.10. 
The relative velocity of collision at the contact is slower than that at the initial time, because of the deceleration due to the 
Coulomb repulsion. 



Motion 


Speed 


Description 


Propagation of charge-equilibrating flow 


0.90 


- 6.5/(0.75x10"^^) fm/s 


Relative velocity for E/A=2.0 MeV 


0.36 


Speed given at the initial time 


Relative velocity for E/A^l.O MeV 


0.23 


Speed given at the initial time 



to lead to fusion or to charge equilibrium for most frag- 
ments. Charge equilibration at higher energies was also 
studied well (for example, see references [8-10, 20-23] of 
Ref. (El). 

As a remark, the previous research on charge equili- 
bration with respect to the collective dynamics is men- 
tioned. As is discussed, charge equilibration was studied 
in association with the isovector giant dipole resonance. 
Charge equilibration is sometimes related to the isovec- 
tor dipole resonance, but not in all cases; i.e. the concept 
of resonance is too restrictive to explain charge equilibra- 
tion. Nevertheless isovector giant dipole resonance means 
that the modes related to the composite nucleus play a 
role, it cannot explain the different upper energy limits 
for different reactions having exactly the same composite 
nucleus. In addition, an isovector mode different from 
the isovector giant dipole resonance sometimes appears 
(see Fig. 4 of Ref. [El). Although it is always true that 
the fast charge equilibration is achieved by the collective 
dynamics, the fast charge equilibration is not necessarily 
achieved by the isovector giant dipole resonance. 

Origin of charge equilibration 

Apart from how and when charge equilibration takes 
place, here we see the reason why charge equilibration 
takes place. First of all, the propagation of zero sound 
is expected to be efficient to any kind of nucleon propa- 
gations and vibrations. The answer is obtained by clar- 
ifying the origin of charge equilibration. When the two 
ions have a contact during the heavy- ion collisions, large 
fluctuations appear in the shape and the internal struc- 
ture. Zero sound is expected to contribute to stabiliza- 
tion by changing both the shape and the internal struc- 
ture, which can be understood by the contribution of each 
term included in the Bethe- Weizsacker mass formula [HI : 

(N — 

B{A) = avol^+asurf^^/^+acoul^^^"^/^+asym^^ 

where the coefficients are given by avoi ^ —16 MeV, 
dsuri ~ 20 MeV, acoui ^ 0.751 MeV and agym ~ 
21.4 MeV, respectively [l6|. After contact, the ffist and 
second terms contribute to stabilization by changing the 
shape such as the volume and surface area, while the 
third and forth terms contribute mainly to stabilization 



by changing the internal structure. It is reasonable that 
the shape change (including vibration) due to zero sound 
propagation leads to stabilization, and we do not go fur- 
ther into detail. Here we focus on stabilization by chang- 
ing the internal structure. The effect due to the symme- 
try energy (the forth term) should play a principal role, 
because the Coulomb energy (the third term) is actually 
small except for the collisions between very heavy ions. 
This symmetry energy is the principal driving force of 
the stabilization by changing the internal structure, and 
its effect acting on zero sound propagation is fast charge 
equilibration. 



SUMMARY 

The propagation of charge equilibrating flow has been 
visualized, and a physical interpretation has been given 
to the fast charge equilibration; i.e., the fast charge equi- 
libration is realized by zero sound propagation. This 
means that the collective dynamics of the many-nucleon 
system and thus its Fermionic statistical properties are 
essential to the fast charge equilibration. In this context 
the upper energy limit for the fast charge equilibration 
corresponds to the limit energy for zero sound propaga- 
tion to be effective. Consequently charge equilibration 
should be regarded as the exchange of nucleons (charge 
exchange), because zero sound does not necessarily entail 
density change. 

As is discussed, fast charge equilibration is not a pro- 
cess appearing only at a certain energy. Fast charge equi- 
libration universally appears in low-energy heavy-ion re- 
actions at energies lower than the upper energy limit, in- 
stead. Here note that the fast charge equilibration pro- 
cess becomes operational only if the two colliding ions 
interact by the nuclear force. 

The correspondence of zero sound propagation in 
femto-scale systems is not only the giant dipole res- 
onance, but also the flow propagating at the Fermi 
velocity. Indeed, the similarity between neutron-rich 
flow shown in Fig. [T] and the isovector giant dipole 
resonance is quite limited. It is suggested that a promi- 
nent role of zero sound in femto-scale systems can be 
experimentally detected by the appearance of the fast 
charge equilibration; whether most of the final products 
of heavy- ion reactions are in charge equilibrium of not. 
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